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Abstract 
At the order of 1/m}, the B meson lifetimes are controlled by the hadronic 
matrix elements of some four-quark operators. The nonfactorizable mag- 
nitudes of these four-quark operator matrix elements are analyzed by 
QCD sum rules in the framework of heavy quark effective theory. The 
vacuum saturation for color-singlet four-quark operators is justified at 
hadronic scale, and the nonfactorizable effect is at a few percent level. 
However for color-octet four-quark operators, the vacuum saturation is 
violated sizably that the nonfactorizable effect cannot be neglected for 
the B meson lifetimes. The implication to the extraction of some of the 
parameters from B decays is discussed. The B meson lifetime ratio is 
predicted as 7(B~)/7(B°) = 1.09 + 0.02. However, the experimental 


result of the lifetime ratio r(A,)/r(B°) still cannot be explained. 
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Heavy hadron lifetimes provide us with testing ground to the Standard Model, 
especially to QCD in some aspects [1-3], because they can be systematically calculated 
within the framework of heavy quark expansion. Theoretically, if we do not assume the 
failure of the local duality assumption, the heavy hadron lifetime differences appear, 
at most, at the order of 1/m¿ [4]. Recent experimental results on the lifetime ratio of 
A, baryon and B meson [5] showed some deviation from the theoretical expectation. 
This has drawn a lot theoretical attentions [6-11]. The current experimental values for 


the lifetime ratios which we are interested in are [5] 


HB) | 

(EN = 1.06 + 0.04 , A 
T(Ap) = 

(BO) = 0.79+0.06. 


This may imply that the Oz) contribution is not enough for the explanation of above 
heavy baryon and heavy meson lifetime difference. To the order of 1/m}, the hadron 
lifetimes have been studied since mid-80s [11, 12, 4, 6, 7]. And the potential importance 
of the 1/m corrections has been pointed out. The parameterization of the hadronic 
matrix elements of four-quark operators which appear in the hadron lifetimes at the 


order of 1/mj is generally expressed as [7] 


(Bloy,(1—15)g7y(1—16)01B) = BiF¿mz, 
(BIb(1 — 75)gg(1 + 75)b] B) = ByFpmE , 0) 
(Blby.(1 — tea’ — ys)tabl|B) = aFam§, , 
(B\b(1 — y5)tagG(1 + Y5)tab| B) = Fam; , 
and m 
Ima load = yag" (1 — ys)blA) = = om , a 
(ABO -aal dl = BEE, 


where the parameters B;, €, (i = 1,2), Fg, r and B should be calculated by some 
nonperturbative QCD method. In above equations, the renormalization scale is arbi- 
trary, and the parameters depend on it. It can be taken naturally at the low hadronic 


scale to apply heavy quark expansion. On the other hand, if the scale is taken at my, 
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parameter Fg(m,) is just the well-defined measurable physical quantity — B meson 


decay constant fp. 


The QCD sum rule [13], which is regarded as a nonperturbative method rooted in 
QCD itself, has been used successfully to calculate the properties of various hadrons. 
In Ref. [8], in the framework of heavy quark effective theory (HQET), the baryonic 
parameters r and B have been calculated by QCD sum rule, r ~ 0.1 — 0.3, B = 1. 
For a complete analysis, the mesonic parameters B; and e; should be also calculated 
from QCD sum rule. The four-quark operators, and hence B;, ci, r and B, are scale- 
dependent quantities when the QCD radiative corrections are included. It was proposed 
by Shifman and Voloshin [11] that at the low hadronic scale, the vacuum saturation 
approximation, namely B; = 1 and e; = 0, makes sense. However in this case, the 
measured lifetime ratio r(A,)/r(B°) cannot be explained [8]. There are some argu- 
ment, on the other hand, that the vacuum saturation maybe a poor approximation [9]. 
Especially from a naive large N. analysis, e;’s are about 1/N. ~ 0.3 [11, 7]. We will 
explore the violation of the vacuum saturation approximation in detail from QCD sum 


rules in the framework of HQET. 


Let us first consider the parameters B;. We construct the following three-point 


Green’s function, 
T° (w, w’) =? I dedye™ ROT (g(x) 7" 5h? OLE, (4 
where w = 2u-k, w! = 2u- k'; h® is the b-quark field in the HQET with velocity v. 


And O denotes the color-singlet operators given in Eq. (2), 
O= brigg ab > (5) 


with DT} = Ta = y*(1 — ys) for Bl and Ti = 1 — y, Ta = 1 + ys for By. In terms of 
the hadronic expression, the parameter B; appears in the ground state contribution of 
10,0), 

Fam 


genen) 


+ resonances , (6) 
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where A = mg — my. The resonance contribution will be simulated by the perturbative 
QCD contribution above some threshold energy due to the local duality assumption. 
On the other hand, this Green’s function r° (w, w’) will be calculated in terms of quarks 
and gluons, that is to say, by the operator product expansion method of QCD. The 
essential feature of the QCD sum rule is that in the QCD calculation, the vacuum 
condensates of quarks and gluons have to be included. Practically the calculation is 
performed at 1 GeV scale or so, only a few terms of the condensates with lowest dimen- 
sions are important. Note that this calculation will be reinforced by the Borel trans- 
formation, and its consistency should be checked through the finding of the so-called 
sum rule window. This procedure will be explained in more detail in the calculation 


of parameters €;. 


The calculation of [°(w,w’) in HQET is straightforward. The fixed point gauge 
[14] is adopted. The assumption that the four-quark condensates are factorizable is 
used. The violation of it will be discussed later. The dominant non-vanishing Feynman 
diagrams are shown in Fig. 1 where the double lines denote the heavy quark. However, 
all these diagrams are factorizable, namely they do not contribute any deviation from 
vacuum saturation. And this is true even when the O(a,) radiative corrections to 
these diagrams are included, simply because O is a color-singlet operator. Note that 
due to the same reason, there is no nonfactorizable gluon condensate contribution. 
The tadpole diagrams in which the light quark lines from the four-quark vertex are 
contracted have been subtracted. Generally, the nonfactorizable diagrams for a color- 
singlet four-quark operator are listed in Fig. 2. They would be the leading contribution 
diagrams which might have an O(1) correction to B; = 1. But it is easy to see that they 
are vanishing because of the special structure of T', and Ts in O (Eq. (5)). Hence, we 
see that the vacuum saturation approximation is valid at hadronic scale for color-singlet 


operators, B; = 1, for i = 1,2, through the leading order consideration. 


To what extent the nonfactorizable effect makes B; deviate from unity? The next 


possibility for non-factorization is to consider two gluon exchanges, like those in Fig. 
3. It is interesting to note that the four-gluon condensate < a?G? > contribution 
to T?(w,w’), like Fig. 3(b), is vanishingf] The perturbative two gluon exchange dia- 
gram, namely Fig. 3(a), can be the leading non-vanishing nonfactorizable contribution. 
(There are some other diagrams which have the same or less order of magnitude numer- 
ically, like Fig. 3(c).) This is a four-loop diagram, its contribution to B; is estimated as 
(248? which is numerically around a few percent and is too small to be important 


for the hadron lifetimes. We obtain that 


B¿=1+0(10%, (i=1, 2). (7) 


The parameters e, (i=1, 2), as we will see in the following analysis, have deviation 
from the expectation of the vacuum saturation at the hadronic scale which cannot be 
neglected for the B meson lifetimes. The procedure is similar to that for parameters 


B;. The three-point Green’s function in this case is constructed to be 
PP (w,w") = 1? f dedye” ANTESALA TO LR, (8) 
where T is the color-octet operators given in Eq. (2), 
T SO! it agql stb. (9) 


In the hadronic language, the parameter e; appears in the ground state contribution of 
DA ua) 
Fame 


ww) = oR eA To) 


+ resonances . (10) 


In the calculation of [7(w,w’), all the condensates with dimensions lower than 6 
are retained. The dominant diagrams are found to be those given in Fig. 4, and 
?This fact is gauge invariant and can be easily observed in the fixed point gauge. In this gauge, the 
light quark propagator with two condensate gluons attached vanishes if one end of the propagator is 


at space-time origin [14]. And in this gauge, the heavy quarks in Figs. 1-3 are free from interaction. 


they are non-vanishing. The four-quark condensate (gq)? diagram vanishes. We have 
neglected the perturbative diagrams which are three-loop diagrams of order a,. That 
means the QCD radiative corrections to e; are not included. From general experience 
of QCD sum rule method, the condensate diagrams play more dominant role than 
the corresponding perturbative one. This neglect is expected viable. The condensates 
parameterize the non-perturbative effects which in the 1 GeV scale are still small 
enough to be treated as power corrections of 1/w and 1/w’ in the operator product 
expansion. Because the perturbative contribution has been neglected, resonances in 
Eq. (10) will be also neglected due to the duality assumption. While the calculation 
can be justified if (-w) and (—w’) are large, however the hadron ground state property 
should be obtained at small (—w) and (—w’). These contradictory requirements can 
be achieved by introducing double Borel transformation for w and w’. There are two 
Borel parameters T and T’ corresponding to w and w’, respectively. They appear 


symmetrically, so we take T = T” in the analysis. 


The sum rules for the parameters e; are 


"i T ei Vf 
a = NG) Zn Fama ’ 
(GG) - T “NT u 
& = -|T + ms (da) Ei 
An 167? Ffm 


There is no gluon condensate contribution to e. Numerically we use the following 


values of the condensates, 
(aa) -(0.23 GeV)? , 


(a,GG) ~ 0.04 GeV‘, (12) 
(970, Gq) = mila), m2~0.8 GeV’. 


K 


For consistence, the HQET sum rule of the parameter Fg [15, 16] will be used. The 
result derived in Ref. [16] is 


Se ie w = 2 
F2mpe "UT = — duv?e-v/* — (qq) h — =] f (13) 
0 


where we is twice of the continuum threshold, which was determined as we ~ 2.0 + 0.3 
GeV. The range of the Borel parameter T in Eq. (11) should be similar to that of 
Fg. This point can be obviously seen from the sum rules for B; if it is written down 
explicitly. In that sum rules, if all the non-factorizable contributions are neglected, 
then the resulting B; is unity only if the Borel parameter is the same as that of Fp. 
Practically we take the window as 0.7 < T < 1.0 GeV. There is no A dependences for 
€; s, because they are canceled actually. The numerical sum rule results for e, and ea 
are given in Fig. 5. From the figures, we see that e;’s have no good stability in the 
window of T. This is because we have not included the perturbative diagram. Finally 


the results for e, and & are obtained as 


¿ 


a = -(41+22)x107?, 
(14) 


K 


6 (6.1 +3.5) x 107? . 


In spite of the large uncertainties in above equation, the numbers are significant for the 
B meson lifetime difference. We see that the vacuum saturation approximation for the 
color-octet matrix elements is indeed violated. The magnitudes of e; can be as large as 


0.1. 


It is necessary to discuss the hypothesis on the four-quark condensate factorization. 
We have used it in above analysis for B; and e;. While this is the working assumption 
for the usual QCD sum rule calculations, the violation of it may imply in our case, new 
contribution to the non-factorizable effect of the four-quark operator matrix elements. 


This contribution to B, and e, is estimated as 


(15) 


where 6(gq)” denotes the deviation of the four-quark condensates from factorization. 
Arguments based on large N, expansion suggest that this approximation is good to 
within 1/N? ~ 10% [17]. Numerically, considering the success of this assumption in 


the calculation of baryons by QCD sum rules [18], we take 9(7q)? to be less than 30% 


T 


of (qq)*. In this case, ôB; and de; are smaller than about 107? which have no significant 


influence on our results for B; and e;. 


It should be noted that the vacuum saturation and its violation we have analyzed 
above are at some hadronic scale, other than the scale of m,. Because we have been 
working in the framework of HQET, in which the natural scale is Uipag K my. The 
renormalization group evolution of the relevant operators was calculated in Refs. [19, 
7]. Information on parameters B; and e; at scale m, is necessary to obtain the hadron 
lifetimes by using the analysis of Ref. [7]. There is a notation difference in this paper. In 
HQET, Fg given in Eq. (13) is in fact a scale-dependent quantity if the renormalization 
effect is considered. From this point of view, Fg in our previous analysis took value 
at hadronic scale, whereas that in Ref. [7] at scale of my. Taking this difference into 
consideration, we have the following evolution relations by choosing as(tinaa) = 0.5 


(corresponding to had ~ 0.67 GeV), 


Bim) ~ Bilttnea) — 0.24€;(Uinad) , 


(16) 
emy) œ —0.05Bi(tnaa) + 072€; (nad) - 
That is 
Bi(m) = 1.01#0.01, Bə(m) = 0.99+ 0.01, 
(17) 
elm) œ —0.08+0.02, em) œ —0.01+0.03, 


from Eq. (14). 


To be more relevant, we discuss a more detailed parameterization for the four- 
quark operators proposed in Ref. [6] in the following. It is motivated for the model- 
independent determination of the hadronic matrix elements from the lepton spectrum 


in the endpoint region of semileptonic B decays. The parameterization is written as 


(Blby, (1 = Y)aayw(l = Y5)b| B) = (Using Vnt = Ising Iv) FM > (18) 


(Blby, (1 u s)taqgy(1 ~~ 15)tab| B) = (VoctUyUy = Jota) Fas > 


where Usingl) Jsingl aNd Voct, Joer are the parameters. The tadpole diagram contributions 
of the matrix elements have been subtracted. The parameters are related to those 


of Eq. (2) as Usingt — 4Gsingt = Bi and Voce — 4goce = €1- In the vacuum saturation 


approximation, Using. = 1, 9singl = Voct = Joct = 0. They can be calculated by QCD 


sum rules. From our above analysis for the color-singlet operator, we have 
Using 1 +0(107%), Ising & O(10°). (19) 


The sum rules for the color-octet operator are obtained as 


(asGG) z le T ¿“NT 
OC = a —— T SD Oo TA 9) 
Voct ( Fr + mo (qq) 8772 Fim, (20) 
E CO nta ) T ne 
— 127 01997) 1672 Fame, ` 
The numerical results for Voce; and Joc: are given in Fig. 6 from which we obtain 
Voce = (1.35 + 0.76) x 107}, gor = (0.74 + 0.40) x 107° . (21) 


They are typically at the order of 0.1. 


From the QCD sum rule point of view, the vacuum saturation for the hadronic 
matrix elements of the color-singlet four-quark operators given in Eq. (5) or in the first 
equation of (18) is valid up to a few percent level. However, that for the color-octet 
four-quark operators is violated at ten percent level. The result that the nonfactorizable 
effect for color-octet operators is more significant can be also understood qualitatively, 
if we look at the perturbative diagrams. This effect is at three-loop level, whereas that 


for color-singlet operators at the four-loop level. 


The 1/m, corrections to above results can be analyzed in principle. While having 
little influence on our above arguments and calculations, they formally belong to the 


O(1/m}) effects to the b-hadron lifetimes. 


In Ref. [6], the extraction of the parameters from B decays are discussed. The role 
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of parameter ggingı might be significant in comparing the decay rates of BP — Iv X,, and 
B- — lvX,,. From our analysis, the difference of the decay rates should be small due 
tO Gsingt is very small even at the scale of my. The ratio of gsingi/ Goce can be obtained 
from the lepton spectrum of the above decays. While having not determined the sign 
of singi, our analysis prefers an even smaller ratio of |Gsingi/Joct| ~ 1/10 than Ref. [6] 


did. These points will be verified by the experiments in the near future. 


The b-hadron lifetime ratios are expressed as [7] 


B- 
¿Ea aa 1+0.03B, — 0.716, + 0.2062 + O(4) , 

TA í (22) 
TA) 098- 0.17e1 + 0.20e2 — (0.013 + 0.022B)r + O(-4) . 

TB") m, 


In above equation, the energy scale for the parameters is at mẹ. the B meson decay 


constant fg is taken as 200 MeV. By using the results in Eq. (17), we obtain 


BT) „1.094002 and TAX 


-CRY (0) ~ 0-98 = 0.01. (23) 


In obtaining above numbers, the QCD sum rule results for the baryonic parameters r 
and B [8] have been also taken values at scale my. From Eq. (23) we see that, while 
T(B-)/T(B?) is agree with measurement, the QCD sum rule result for (A,)/r(B°), 


after including 1/mj corrections, still contradicts the experimental measurement. 


In summary, the nonfactorizable contributions of the hadronic matrix elements of 
four-quark operators relevant to the B meson lifetimes have been studied by QCD 
sum rules in the framework of HQET. The vacuum saturation for color-singlet four- 
quark operators is justified at hadronic scale, and the nonfactorizable effect is at a few 
percent level. However, the vacuum saturation for color-octet four-quark operators is 
violated sizably that the nonfactorizable effect cannot be neglected for the B meson 
lifetimes. The implication to the extraction of some of the parameters from B decays 
has been discussed. The B meson lifetime ratio has been predicted to be consistent 


with experiment, 7(B~)/7(B°) = 1.09+0.02. However, the discrepancy between theory 
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and experiment on r(A,)/r(B°) has not been solved. It is unlikely that higher 1/m» 
corrections will give the solution. If the experimental data of r(A,)/r(B°) is further 
confirmed in the future, that may imply either the failure of local duality assumption 


or some new physics. 
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Figure captions 


Fig. 
Fig. 
Fig. 
Fig. 
Fig. 
Fig. 


. Dominant non-vanishing Feynman diagrams for T? (w, w”). 

. Nonfactorizable diagrams for a general color-singlet four-quark operator. 

. Nonfactorizable diagrams by two gluon exchanges for T° (w, w’). 

. Condensate contribution to [7 (w, w’). 

. Sum rules for e, (a) and ez (b). The sum rule window is T = 0.7 — 1.0 GeV. 


. Sum rules for Vo. (a) and goet (b). 
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Figures 


Fig. 2 


(b) 


Fig. 4 


Fig. 5(b) 
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Fig. 6(b) 
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